Velocity, temperature and normal force dependence on friction: 
An analytical and molecular dynamic study 
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In this work we propose an extension to the analytical one-dimensional model proposed by E. 
, Gnecco (Phys. Rev. Lett. 84:1172) to describe friction. Our model includes normal forces and 

/"*\ ■ the dependence with the angular direction of movement in which the object is dragged over a 

surface. The presence of the normal force in the model allow us to define judiciously the friction 
C^") , coefficient, instead of introducing it as an a posteriori concept. We compare the analytical results 

with molecular dynamics simulations. The simulated model corresponds to a tip sliding over a 
surface. The tip is simulated as a single particle interacting with a surface through a Lennard- Jones 
(6 — 12) potential. The surface is considered as consisting of a regular BCC(OOl) arrangement of 
particles interacting with each other through a Lennard- Jones (6— 12) potential. We investigate the 
f) |' system under several conditions of velocity, temperature and normal forces. Our analytical results 

are in very good agreement with those obtained by the simulations and with experimental results 
from E. Riedo (Phys. Rev. Lett. 91:084502) and Eui-Sung Yoon (Wear 259:1424-1431) as well. 



PACS numbers: 46.55.+d, 07.79.Lh, 07.79. Sp, 81.40.Pq, 62.20.Qp, 68.35. Af 
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^ ! I. INTRODUCTION 

Understanding the origin of tribological phenomena is a fascinating and challenging cnterprize. The classical point 
of view of the frictional phenomena, can be synthesized in the three laws of friction, valid in the macroscopic scale 

' 1. Friction is independent of the apparent area of contact, 

2. Friction is proportional to the applied load. The ratio between the friction force and the applied load is named 
£Nj ■ the coefficient of friction (/i = /l//jv) and it is larger for static friction than for kinetic friction, 

I | 

[ — ' 3. Kinetic friction is independent of the relative sliding velocity. 

O ' 

Since new tools, such as the atomic force microscopy (AFM), have made possible to examine the friction phenomenon 
in great detail these laws have been questioned in systems with dimensions approaching the nanometer scale. At 
the same time, the development of ultra fast computers have allowed to test new theories on the nano-scale friction 
5^ ' world. Although tribology is an old science, and in spite of the efforts and progress made by scientists and engineers 
in the last years, tribology is still far from being a well-understood subject. In fact, it is incredible that even knowing 
several properties as surface energy, elastic properties and loss properties, a friction coefficient cannot be found by 
using an a priori calculation. Although in the macroscopic scale the friction force, /l, is independent of the relative 
velocity, in the nanometric scale some authors [H, 0, HI, @] observed that the mean value of the friction force presents 
a logarithmic velocity dependence. Another important result was the conclusion that friction force is proportional to 
the effective contact area down to the nanometer scale 0j- An analytical one-dimensional model known as Tomlison 
model [l[ was able to explain several features of the nanoscopic friction. Using the Tomlinson model in the limit of 
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low velocities Gneco et al Q showed that the friction force has a logarithmic dependence with the velocity. Using 
the same ideas, but in the limit of higher velocities, Sang et al. Q obtained that the friction force is proportional 
to | ln(v)| 2 / 3 , were v is the relative velocity. Using a first principle model Persson Q was able to show that in the 
limit of small contact areas the result of Sang is recovered while in the limit of large contact areas the Gneco result 
fitted better. The aim in this work is to develop a model from first principle by extending the one-dimensional model 
proposed by Sang et al. [H to three-dimensions. Based in our approach we obtain a friction coefficient which can 
be calculated knowing simple parameters of the model (As bound energies and the positions of the minima between 
atoms.). Such parameters can be obtained by using ab-initio calculations or measured experimentally by using FFM 
(Friction force microscopy) [1, [H, [l(| 03 or DFS (Dynamic force microscopy) [l2[ . We study the sliding frictional 
process by using two approaches: 

• Developing an analytic model that considers the potential energy between a atom in the tip and the surface 
atoms described by the model presented by W. A. Steele [l3| as a sum of pair-wise (6 — 12) LJ potentials. The 
analytic treatment extends the one dimensional model proposed by Riedo et al.0] including the normal force 
and as a consequence, the effects of adhesion energies. 

• Using MD simulations by considering that the potential energy between the tip's atom and the surface atoms 
is described as a sum of (6 — 12) L J potentials (See ref. [l4| and references there in.). 



II. FIRST PRINCIPLE MODEL 




Pi gure 1: Perspective view of a typical initial configuration [BCC(OOl) geometry] for the FFM experiment. 



In this section we show a general picture of the nature of the kinetic friction. Consider the sliding system shown 
in figure |T]) in which one particle of mass m is connected through a spring to a cantilever or drive. The particle 
experiences a total force described by the potential [l[ 
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represents the harmonic spring constant complying the cantilever with the tip, Vi„t(q) is the surface-tip corrugated 
potential, q = {q x ,<ly,<lz) are the coordinates of the tip and r = (x,y,z) the coordinates of the support. As we are 
interested to study the influence of the normal force on the sliding process, let us first note that the critical state 
(Where we denote the drive position by r c = (x c , y c , z c ) and the particle position by q c = {q xc , qyc, Qzc)) is the position 
where the tip jumps from a stable position in the surface to the next one. To illustrate the occurrence of the critical 
state, we show in the figure [21 the total potential energy of the tip as a function of q x , for two different positions of the 
cantilever, x < x c (full line) and x = x c (dashed line). As a matter of clarity we restrict this figure to the x direction. 
From figure [H the critical point is defined as the inflexion point of the total potential energy and mathematically it 



means, 
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Figure 2: Schematic figure illustrating the critical state for two different positions of the cantilever, x < x c (full line) and x = x c 
(dashed line). The marked points define the energy barrier that prevents a jump from one stable position (q x -) to another 
(equation [8}. Figure from reference Q- 
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The first condition, equation ([3]), is always satisfied at equilibrium, it states that the total force on the particle must 
vanish. The second condition, equation is satisfied in a transition from a stable to an instable position on the 
plane of the surface (critical points). It states that the determinant of the Hessian matrix in the plane vanishes. It 
follows from the fact that at these points, the slope of the force due to the substrate equals the slope, k a , of the spring 
force. The third condition, equation ([5]), defines that the particle is always in contact with the surface. 

We can write the total potential as Vr t(<f, t) = Vr t(qj| , r]| , q 2 , z)i where qu and rj| are coordinates parallel to the 
surface and the others correspond to the normal components. We can expand the potential around the critical points, 
X = {Qc,r c ), as 
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As a matter of clarity we write this potential in a more compact form, 

V To t ~ A(x, y, z) + Ci{x - x c )(q x - q xc ) 

+ 2 (y- y c )(q y - q yc ) + C 3 (z- z c )(q z - q zc ) 

+ Ci 3 (q x - q xc ){q z - q zc ) + {E 3 /2)(q z - q zc ) 2 
+ C 23 {q y - q V c)(qz - qzc) 

+ (E 1 /6)(q x - q xc f + {E 2 /Q)(qy ~ q yc f + ■■■ 



(6) 



(7) 



where the constants, A, C\, C 2 , C 3 , C13, C 23 , Ei, E 2 and E 3 will be explicitly written in the next section, they 
are immediately obtained by direct comparison with equation The main goal of this model is to show that the 
dynamic energy barrier, AE, which prevents the jump of the tip from a stable equilibrium position on the surface to 
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the nearest neighboring one, is a function of the instantaneous lateral force, /l, and the instantaneous normal force, 
f z . The energy barrier AE is defined as 



AE — VrotHnax, t) — VTotyQmini t) 



(8) 



The quantities q max and q m i n correspond to the first minimum and maximum of the combined potential at time, t, 
it is determined by the equilibrium condition dV/dq = [J] (See figure [5]) . 
Using the equilibrium condition in the equation J7J we obtain, 



E 



Jrfax ~ Ixc) 2 + Ci 3 (q z - q zc ) + Ci(x - x c ) = 
-7r{% ~ lycf + C 23 (q z - q zc ) + C 2 (y - y c ) = 



(9) 
(10) 



E 3 (q z ~ qzc) + Ciz{q x - q xc ) + C 23 (q y - q yc ) 

+C 3 {z-z c ) = (11) 

substituting (fTTj) in §§§ , (fit))) and holding only terms of first order in C\ 3 and C23 , we obtain two quadratic equations 
in q x , q y to solve, whose solutions are 



with 



q x ± ~ q xc = ±—^^2C 1 {x' c - x) 
E i 

q y ± - q V c = ±-^^2C 2 (y' c -y) 

E 2 

dx c = C lXc - {C 13 /E 3 )C 3 (z c - z) 
C 2 y'c = C 2 y c - (C 23 /E 3 )C 3 (z c - z) 



(12) 
(13) 



Using the definition of AE, and equations from ([7|) to (fT3|) we can get the energy barrier as 

2l\[C 1 {x c -x)fl* , [C 2 (y' c - y)] 3 / 2 1 



AE = 



1/2 



E. 



1/2 



Defining, 



Ut) = C\x(t) 
f y (t) = C 2 y(t) 
fS) = C 3 z{t) 



^xc C\X C , F xc — C\X C 
Fyc = C 2 y c ; F yc = C 2 y c 
F zc — C 3 z c , 



/*(*) = /i(*)cos(0) ■ f y (t) = f L (t)sm(6) , 
F xc = F lc cos(6) ■ F yc = F lc sm(9) , 

Flo = s/Fl + Fj c ; f L (t) = yff x (t) 2 + fy(t)* 
f' xc = F; c cos(6) ; F yc = Fism(6) , 
we find an energy barrier which is similar to that found in other works [H, [^, [TlT ] , 

3/2 
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FiAfS)) = 



Mi* + - Fzc) 2 + 2F lc (n 0x cos(0) + sin(0))(/ 2 (i) - F zc ) 



(18) 



We note that AE presents a dependence with the instantaneous normal force, f z , and the angle < 6 < ir between 
the direction x and the force pushing the cantilever. Now we will study the consequences of this dependence in 
the friction force. At zero temperature, the lateral force required to the tip jumping from a minimum to another is 
/l(£) = F lc (f z (t)) (i.e., the force corresponding to AE = 0). At finite temperature, T, the occurrence of thermally 
activated transitions between two minima, when AE > 0, leads t o fri t) < F lc (f z (t)) and the probability that the tip 
does not jump, p(t), is described by the master equation [1,0,11, 
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(19) 



where [3 = 1/kbT, isq is the jump frequency transition and kb is the Boltzmann constant. Using the condition for the 
maximum jumping probability, d 2 p/dt 2 = 0, the following expression has to be satisfied: 
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and K^ff is the effective stiffness of contact [15| in the directions a = x,y, z. It is important to stand out that f\,ff 
is a function of the normal force. This dependence, will appear in a crucial way in the deduction of the theoretical 
coefficient of friction that we will present. The effective stiffness of contact can be experimentally determined by 
measuring the inclination of the "stick region" , which is the region were the tip is in a stable position before the 
equation (0| is satisfied in the friction force loop [l6j]. Thus, the effective stiffness of contact is defined as, 
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where K^ ontact is the stiffness of contact which is a function of the normal force [13, 13 and k x (equation @), is the 
spring constant of the cantilever. Using equations (fTS)) . (|2H)) and (12"T]) we obtain the lateral force, 
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This result has a similar form to those presented by Riedo et. al 0] with three interesting physical differences: 1) 
the explicitly normal force dependence on F [c (f z ) and t>o(/z) as indicated by experimental results 0]; 2)the angular 
dependence of the friction force relative to the periodic crystalline structure of the sliding surface; 3) the logarithmic 
normal velocity dependence of the friction force, as Jeon et. al. [l9j have shown experimentally and numerically 
simulated using a single particle model. Note that by introducing a normal force dependence in the friction force it 
is natural to think of the coefficient of friction as defined by 
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where 



dip dn x eff 
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Observe that this equation gives a friction coefficient that is a function of velocity, temperature, normal force and 
the direction that the cantilever is dragged in the surface. In the next sections we show this results for a selected 
potential and we compare it with MD simulations. It is important to note that if we know the potential between the 
surface and the tip atoms we can predict the value of the friction coefficient with no adjustable parameters. 

In the next section we will use a total potential to calculate the values of the parameters defined in the model and 
show that the theoretical friction force agree very well with experimental results [3, 0] ■ 



III. SURFACE-TIP'S POTENTIAL AND PARAMETERS OF THE ANALYTIC MODEL 



The model we developed has no adjustable parameters. Once the potential energy is given, we are able to calculate 
the friction force. For a well behaved surface, Steele[l3| has derived a model potential which is very fair in describing 
the interaction with a periodic surface 

V int (q\i,q z ) = V (q z ) + Fife) ^ cos((? • q {l ), (30) 

G 

where q\\ = {q x ,qy) are the coordinates of the tip's atom parallel to the substrate and G is the set of the six shortest 
reciprocal lattice vectors of the substrate. The first term in equation ([30)) describes the mean interaction of the atoms 
with the substrate, and the second term describes the periodic corrugation potential. Expressions for Vo and V\ 
were derived by Steele [l^], assuming that the substrate potential Vi n t(q) is a sum of LJ potentials between one film 
atom and all of the atoms in the substrate. The parameters of the LJ potential can nowadays be experimentally 
determinated by DFS [2(| experiments for a specific material. Other expressions for Vo and Vi, where described by 
some authors such as Persson et. al. [2l|, Liebsch ct. al. [13] that consider these terms composed by exponentials 
and Tomassone et. al. p3| consider expressions that give a correct description of the interaction of a metallic surface 
with a noble gas atom. In this work we use the expression derived by Steele. At this point we need to define the 
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surface arrangement of particles which we are interested to work with. We use the regular BCC(OOl) surface, so that 
the total potential energy (Equation (Q])) becomes 
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Solving these equations we get, 
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As we can see, equations Q36p and (|37p define a transition from a stable to instable tip's position, or in other words 
a transition from the stick to the slip state. Using this information we are able to explicitly calculate the parameters 
defining the friction force of the analytical model presented above. With this in mind, we explicitly define and calculate 
the parameter C\, C 2 , C 3 , Ci 3 , C 2 3, Ex, E 2 and E 3 which led to equation (JT]). From equation (|3"Tj) we obtain, 
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Using equations ([32 ]) . (|33 ]l . ([3i ]) . (|36 |) and t[37 j) we get 
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To calculate \i$ and we start by calculating C13, C23 and E 3 . From equations (f3"Tj) . and (j3"7) we can 

write 
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Using these results in equations (jT5J) we obtain, 
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From equations (f3"Tj) . (|36|) and (l37|) we obtain id and E 2 by, 
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where as defined by equation ([T7]) we can calculate Xi(0) as, 
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Note that all the parameters Fi c , F* , fiQ, and Ai are functions of the critical point q zc . Using the functions Vo(q z ) 
and V\(q z ) of the Steele potential defined by fl3| . 
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where p is an integer, q is the total number of atoms per unit surface cell, Aq z is the distance between planes, a s = af, 
is the area of the unit lattice cell, a\ = <t ss , a ss is the nearest neighbor in the solid and at s and £t s are the parameters 
of a (6 — 12) Lennard-Jones potential between the tip and the surface atoms, and 
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where K n is the modified Bessel function of second kind and g\ = G x = G y = 2tt ja\. Using the equations from (f50|) 
to ([55]) we plot in figure [3] the parameters fi 0x = /j, , and Ai and Vi nt (q» = (ai,a\),q z ) as a function of q z , for values 
of rj x = rj y < 1 since we are interested in real values of the parameters. 



A. Theoretical prediction of the friction force and coefficient of friction 



In the figure [3] we show the parameters fj,Q X , X% and Vi nt calculated by using the equations (|4"5)) . (|4"9")l and (|3"0|) 
respectively. Using those results we can calculate the friction force, (/l), and the coefficient of friction, (/i). The 
figure d] shows (/^) and (/i) as a function of the relative velocity for several values of the of the normal force, (/ 2 ). 
We note that these results agree quite well with the experimental data of Riedo et.al 
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Figure 3: The parameters fiox = Mo, Ai and Vi„t(q\\ = (ai,ai),q z ) plotted as a function of q~. Here 9 — 0, a ts = 1.2a as , 
sts = 0.5£ ss (full), 0. le ss (dashed). 




V L 



Figure 4: Theoretical friction force as a function of the velocity for different values of the normal force given by the equation 
(|25|) . Here from bottom to top are shown /jvi < /at2 < /at3 < /at4 < /ats respectively. 



In the figure [5] we present the dependence of the friction force (left) and the coefficient of friction (right) with the 
normal force (/ z ). We note that the dependence of the coefficient of friction (right) with the normal force (f z ) agree 
quite well with the experimental data of Eui-Sung Yoon et. al Q . 

In the figure [6] we present the dependence of the friction force (left) and the coefficient of friction (right) with the 
pushing angle, 9, relative to the direction (100) of the lattice. It is important to note that this result is a particular 
case for the potential defined by the equation (|3l"j) . It is expected that the effect strongly depends on the symmetries 
of the surface. 



11 




Figure 5: Friction force(left) and friction coefficient(right) as a function of the normal force for different values of the velocity. 
Here from bottom to top are shown Vl = 0.045, vl = 0.25, vl = 0.5, vl = 0.75. 
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Figure 6: Dependence of the friction force(left) and the friction coefficient(right) on the direction of the sweepings, 6, relative 
to the (100) direction on the surface. 



IV. SIMULATION 



In this section we present our simulation for the system discussed above. The simulation is carried out using 
molecular dynamics (MD). The simulational arrangement is as shown in figure [TJ which represents a tip interacting 
with a surface. The tip is simulated by a single particle. Three springs are attached to the particle, two parallels 
to the surface plan (x,y directions) and the other perpendicular to the surface plan (z direction). This arrangement 
allow us to measure the forces, normal (f z ) and parallel to the surface (f x ,fy)- The surface is represented by an 
arrangement of particles which interact with each other through a Lennard-Jones (LJ) (6 — 12) potential 



J V hJ > 



fajiuj) - <Pi,j(r c ) - (n.j - r c ) r'",,,." ) if nj < r c 

\ z '3> / nj—r c 

if r hj > r c 



(56) 



/ \ 12 




fe) - 









The indexes i and j stands for particles at position n and rj respectively, and 1 < i,j < N, where is the total 
number of particles, Oi^ = a ss , = £ss between surface atoms and ai.j = a ts a ss , tij = e^ess between the tip and 
surface atoms. A cutoff is introduced in the potential in order to accelerate the simulation. If the force on a particle is 
found by summing contributions from all particles acting upon it, then this truncation limits the computation time to 
an amount proportional to the total number of particles A^. Of course, this truncation introduces discontinuities both 
in the potential and the force. To smooth these discontinuities we introduce the constant term 4>(r c ). Another term 
{d(j>i t j(rij)/drij) r _ r is introduced to remove the force discontinuity. Particles in the simulation move according 
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Newton's law of motion, which generates a set of 3N coupled equations of motion which arc solve d by increasing 
forward in time the physical state of the system in small time steps of size St = 10~ 3 cr ss yj (m/e ss ). The resulting 
equations are solved by using Bccman's method of integration. In order to improve the simulations we use a Vcrlet 
table and cell division method. The temperature, T, of the surface can be controlled by using a velocity renormalization 
scheme [13, US US H3, HE HI, US HH . From the equipartition theorem we can write that 



3^T. 

m 



(58) 



We want to control the value of (v 2 ) to correspond to a chosen temperature Tf. By initializing the system with (i |2 )c 
we multiply each velocity by a factor ao 



a = 



3k B T 



f 



(59) 



By evolving in time the system we can create a sequence {a n }, such that after a finite number of time steps the 
temperature of the system converges to Tf. We measure the time t, and temperature T, in units of a ss y/mj 'e ss and 
^ss/ks respectively. 



A. Numerical background 



Our simulation is as follows. We consider the system as consisting of an arrangement of particles of mass m, 
coupled by the Lennard- Jones potential defined by the equation (|56|) . The system is arranged in 4 layers with free 
boundary conditions in all directions. The first layer is frozen in a regular arrangement as in the (001) surface of 
a Lennard- Jones bec crystal in order to maintain the whole structure as flat as possible. The tip is simulated as a 
single particle of mass m, attached to three springs of elastic constant k x — e ss /a 2 s ,k y = e ss /a 2 s and k z = e ss /a 2 s 
as shown schematically in figure Q] With the tip close to the surface we thermalize the system at temperature T. 
After thcrmalization the tip is pushed in a direction parallel to the surface at constant velocity vq x = «lcos(#), 
vq v = vl sin(#), vqz — 0, vl = v^J (e ss /m). Here 9 is defined relatively to the x direction and v is changed to obtain 
different velocities. For each simulation the distance between the tip and the surface is fixed at the beginning of the 
process, so that we can control the perpendicular force on the tip. By measuring the size variation of the springs we 
can calculate the laterals, f x , f y forces and the perpendicular, f z force on the tip. This forces are measured in units 
of, e ss /a ss . The velocity, position, energy and forces are stored at each time step for further analysis. Before we start 
the simulation wc have to estimate the melting temperature T m of the system. Figure [7] shows the total energy (per 
particle) as a function of temperature. The melting temperature is estimated as the inflection point of the curve. 
We find T rn f=a 1.1 in accordance with earlier calculations [13, [H, HfJ- Based in this result our simulations will be 
performed with T < T m and the temperature will be specified in each results. We have simulated the system for 




Figure 7: 
T m w 1.1 



Energy as a function of temperature. The melting temperature is estimated as the inflexion point being around 



several velocities, temperatures, initial distance of the tip to the substrate or equivalently, normal force in the tip and 
dragging angles relative to the surface. We are mainly interested in studding the effects of the velocity, the normal 
force and dragging angles relative to the surface in the friction force and coefficient of friction. 
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B. Dependence with velocity and normal force 



In this section we show the results of our simulations for several values of the relative velocity tip-surface and for 
five different values of the normal force. In all the MD simulations of this subsection we use = 0, T = 0.5. In the 
figure |(5J) we show the simulation (points) and the theoretical (full lines) results for the friction force as a function 
of the velocity for some normal forces. Note that our theoretical result, equation ([25]) . is in very good agreement 
with the results of the MD simulations, only adjusting appropriate values for the critical position of jumps, q zc (See 
section HU . In the figure (J9j) we show the results of the simulation (points) and the theoretical results (full line) for 
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Figure 8: Friction force as a function of the velocity for some selected values of the normal force. From a) to e) we have 
respectively < f, >~ 0.95, 1.54, 2.25, 2.75 and 2.50. The points are the results of the simulations and the lines are from our 
analytical approach. 
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Figure 9: Averaged friction coefficient as a function of the velocity, calculated from the figure [8] The points are the results of 
the simulations and the line is the analytical result. 

the averaged friction coefficient as function of the velocity. These results are in good agreement with the simulation, 
indicating that our model is consistent with the simulation. 



14 



2 
1.8 
1.6 

f 14 
Jl 1.2 

1 

0.8 
0.6 



- 


~ 


_ 












i ' iii 


i ■ n 


0.1 


1 












Figure 10: Friction force as a function of the velocity for some selected temperatures. The results of the simulations are shown 
as circles the lines are from our analytical approach. The results from a) to f) are respectively to T = 0.1, 0.2, 0.3, 0.4, 0.5 and 
0.6. 




Figure 11: Averaged friction coefficient as a function of the velocity for some selected temperatures. Here, the lines are for the 
analytical results and from bottom to top corresponds respectively to T = 0.1, T = 0.2, T = 0.3, T — 0.4, T — 0.5 and T — 0.6. 



C. Dependence with Temperature and Velocity. 

In this section we present some simulation for the friction force for several velocities, and different values of tem- 
perature, (T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6). In all the simulations in this subsection we use 9 = and < f z >~ 0.95. 
In the figure (fT0|) are the results of the MD simulations (points) and the theoretical result (full lines) for the friction 
force. Note that our theoretical results arc in excellent agreement with the results of the MD simulation. In the figure 
(fTTjl we show the analytical results for the averaged friction coefficient as a function of the velocity. In this figure, the 
lines from bottom to top corresponds respectively to T = 0.1, T = 0.2, T = 0.3, T = 0A,T = 0.5, T = 0.6. 

In figures (J9j) and (|11| , it is shown the behavior of the friction coefficient as a function of the relative velocity. We 
note that initially the friction coefficient increases with the velocity until it reaches a maximum and then starts to 
diminish. The initial increasing and the maximum value for the coefficient of friction can be related to the fact that 
when the velocity increases the tip executes jumps more and more energetic, and a high friction mode associated with 
single slips appears as indicated by Nakamura et. al. [32$ . As the velocity pass though a certain limit a low friction 
mode with double slips appears leading to the observed reduction of the friction coefficient. 
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Figure 12: Friction force as function of the angle, 6, for some selected normal forces. From a) to e) we have respectively 
< fz >~ 0.95, 1.54, 2.25, 2.75, 2.50. The points are the results of the simulations and the lines are our theoretical results. 

D. Dependence with Normal Force and Force Pushing Angle. 

Below we present the simulations for five different values of the normal force and several pulling angles, < 9 < ir/2, 
between the direction x and the force pulling the cantilever. In all the MD simulations in this section we use vi = 0.05, 
T = 0.5. In the figure (fP2")l and (fT3)) we show our results of the simulation (points) and the analytical calculations 



^0.25 




0.2 



Figure 13: Averaged friction coefficient as a function of the angle, 6, calculated from the figure [121 The points are the results 
of the simulations and the line is the analytical result. 



(full lines) for the friction force and the averaged friction coefficient as a function of the angle, 9. Note that our 
theoretical result, equation (|25|) . is in very good agreement with the results of the MD simulations. 
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V. CONCLUSION AND COMMENTS 



In this paper we have proposed a tri-dimensional model to describe friction force. Our approach generalize the 
one-dimensional friction force model proposed in earlier works 0, H[ . In our model we include the influence of the 
normal force which allow us to calculate the friction coefficient for first principles. We apply our model to calculate 
the friction force and friction coefficient for a tip-surface interaction. The interaction between the tip and the surface 
is represented by a Lennard- Jones potential and the surface represented by a potential that represent a (001) surface 
symmetry and the friction force obtained are in good agreement with the experimental results presentes by RiedoQ. 
With the intend to test our model we performed classical molecular-dynamics simulation for a particle attached in 
a spring and pushed over a BCC(OOl) surface crystal. For this purpose we performed the simulation for several 
velocities, temperatures, normal forces and dragging angle. We have found that the friction forces and the friction 
coefficient are in good agreement with our molecular dynamic results, proving the effectiveness of our model. 
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